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Abstract 

We give a new characterization for the convergence in distribution to 
a standard normal law of a sequence of multiple stochastic integrals of 
a fixed order with variance one, in terms of the Malliavin derivatives of 
the sequence. We also give a new proof of the main theorem in [7] using 
techniques of Malliavin calculus. Finally, we extend our result to the 
multidimensional case and prove a weak convergence result for a sequence 
of square integrable random variables. 
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1 Introduction 

Consider a sequence of random variables Fk belonging to the nth Wiener chaos, 
n > 2, and with unit variance. In [7], Nualart and Peccati have proved that this 
sequence converges in distribution to a normal N (0, 1) law if and only if one of 
the following two equivalent conditions hold: 

i) limfe^oo E(F4) = 3, 

ii) limfc^oo fk iXi; fk = 0, for all 1 < Z < n - 1, 

where fk is the square integrable kernel associated with the random variable 
Fk, and fk ®; fk denotes the contraction of / indices of both kernels. In a 
subsequent paper, Peccati and Tudor [lOj gave a multidimensional version of 
this characterization. 
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There have been different extensions and apphcations of these results. In [3] 
Hu and Nualart have appHed this characterization to estabhsh the weak conver- 
gence of the renormahzed self-intersection local time of a fractional Brownian 
motion. In two recent papers, Peccati and Taqqu [HI study the stable con- 
vergence of multiple stochastic integrals to a mixture of normal distributions. 

The aim of this paper is to provide an additional necessary and sufficient 
condition for the convergence of the sequence Fk to a normal distribution, in 
terms of the derivative of in the sense of Malliavin calculus. This new 
condition is 

WDFuWh , -> n. (1) 

On the other hand, we give a simple proof of the fact that condition ([T]) implies 
the converge in distribution to the normal law based on Malliavin calculus. The 
main ingredient of the proof is the identity SD = — L, where 5, D and L are the 
basic operators in Malliavin calculus. In this way, we are able to show the result 
of Nualart and Peccati with the additional equivalent hypotheses ([Ij , without 
using the Dambis-Dubins-Schwartz characterization of continuous martingales 
as a Brownian motion with a time change. 

We also discuss the extension of these results to the multidimensional case. 
In the last section we study the weak convergence of a sequence of centered 
square integrable random variables. The result assume conditions on the Malli- 
avin derivatives of the chaotic projections of the sequence {Ffej^gN ^^"^ 
derivatives of the sequence . Therefore, it can be used with non regular 

random variables in the Malliavin sense. 

In [8] Peccati and Taqqu also provide a condition for convergence to the nor- 
mal law, involving projections of Malliavin derivatives. This type of condition 
is different from ours, and it is inspired on Clark-Ocone's formula. 

Condition ([T]) is a useful tool in establishing the central limit theorem for 
sequences of random variables defined in terms of a fixed function of a Gaussian 
process. We apply this approach to derive the weak convergence of the normal- 
ized sums of odd powers of the increments of a fractional Brownian motion. In 
this case, condition ([1]) follows easily from the ergodic theorem. 

The paper is organized as follows. In Section 2 we introduce some notation 
and preliminary results. In section 3 we state and prove the main result of the 
paper. Section 4 deals with the multidimensional version of the result proved in 
Section 3, and in Section 5 we apply the previous results to the weak convergence 
of a sequence of centered and square integrable random variables. Finally, in 
Section 6 we discuss the application of our approach to an example related to 
the fractional Brownian motion. 

2 Preliminaries and notation 

Let H he a, separable Hilbert space. For every n > 1 let H^" be the nth tensor 
product of H and denote by ff®" the nth symmetric tensor product of H, en- 
dowed with the modified norm -v/n! IMIji/on. Suppose that X — {X (h) : h ^ H} 
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is an isonormal Gaussian process on H. This means that X is a centered Gaus- 
sian family of random variables indexed by the elements of iJ, defined on some 
probability space (f2,^, P), and such that, for every h,g £ H, 

E[X {h)X{g)]^(h,g)H. 

We will assume that is generated by X. 

For every n> 1, le Hn be the nth Wiener chaos of X, that is, the closed linear 
subspace of (O, T, P) generated by the random variables {Hn {X (h)) ,h G 
^A\^\\h ~ 1}' where i?„ is the nth Hermite polynomial. We denote by Ho 
the space of constant random variables. For n > 1, the mapping In{h'^'^) = 
n\Hn {X (h)) provides a linear isometry between if®" and Tin- For n = 0, 
Hq — M, and Iq is the identity map. 

It is well known (Wiener chaos expansion) that {il, !F, P) can be decom- 
posed into the infinite orthogonal sum of the spaces 'Hn- Therefore, any square 
integrable random variable F G {^,^, P) has the following expansion 

oo 
n=0 

where /o = E[-F'], and the /„ G 7J©" are uniquely determined by F. For every 
n > we denote by J„ the orthogonal projection on the nth Wiener chaos Hn, 

so In if n)=Jn{F). 

Let {cfc, fc > 1} be a complete orthonormal system in H . Given / e -ff®"and 
g S i?®"*, for Z = 0, ...,nAm the contraction of / and 5 of order / is the element 
of iJ»("+™-2;) defined by 

00 

f®ig= ^ (/, (g) • • • (g) eijj:^g,i (g) (g, (g) • • • (g) ei,)^®i . 

ii,...,ii=l 

We denote by f®ig its symmetrization. Then, f ®o g = f ^ g equals to the 
tensor product of / and and for n ^ m, f (g„ g — (/, g)^®^. 

Let us introduce some basic facts on the Malliavin calculus with respect 
the Gaussian process X. We refer the reader to Nualart [H for a complete 
presentation of these notions. Consider the set S of smooth random variables 
S of the form 

F^ f{X{h,),...,Xihn)), (2) 

where hi, hn £ H , f G (R") (the space of bounded functions which have 
bounded derivatives of all orders) and n G N. The derivative operator D on a. 
smooth random variable of the form ([2|) is defined by 

DF = J2§^iX{hi),...,X{hn))h,, 

i—l * 

which is an element of (fl; H). By iteration one can define D"''F which is an 
element of (il;i7®™). For m > 1 we denote by D™'^ the completion of S 
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with respect to the norm ||^||^ 2 E^^^^ by 

1=1 

We denote by S the adjoint of the operator D. That is, 5 is an unbounded 
operator on {fl; H) with values in (fi), whose domain, denoted by Dom(5, 
is the set of i7-valued square integrable random variables u € LP' {VL\ H) such 
that 

n{DF,u)H]\<c\\F\\^,^^^, 

for all F E D^'^. If u belongs to Dom(5, then S {u) is the element of {fl) 
characterized by 

E[Fd iu)]=n(DF,u)H], 

for any F G D^'^. 

The operator L defined on the Wiener chaos expansion as L = X^i^o 
is called the infinitesimal generator of the Ornstein-Uhlenbeck semigroup. The 
domain of this operator is the set 

oo 

DomL = {FeL^ {ft) : \\JnF\\l2^n^ < +00} = D^-^. 

n=l 

The next proposition explains the relationship between the operators D, 5 and 
L. 

Proposition 1 For F E (il) the statement F E DomL is equivalent to 
F E DomSD (i.e., F E D^'^ and DF E DomS), and in this case 

5DF = -LF. 

In the particular case where H — {A, A, /i), {A, A) is a measurable space, 
and /J, is a fi-finite and non-atomic measure, then H®" = L'^ [A"' , A^"' , fj.^^) 
is the space of symmetric and square integrable functions on A" and for every 
/ E jff®", /„ (/) is the multiple Wiener-Ito integral (of order n) of / with respect 
to X, as defined by ltd in In this case, F — '}2^=a ^nifn) & ^^'^ if and only 
if 

00 

n=l 

and its derivative can be identified as the element of L^{A x il) given by 

C30 

Af^ = ^«/„-i(/n(-,i)). (3) 

n=l 

We need the following technical lemma. 
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Lemma 2 Consider two random variables F ~ Inif), G = Im{g), where 
n,m > 1. Then 



E 



nAm , I |n2 

{DF,DG)l] = J, ,u \u Un^ ll/^^5||Lo,..™-.., ■ (4) 

{{n - ry.[m ~ ry. [r - ly.) 



Proof. Without loss of generality we can assume that H = {A^ A, /i), [A, A) 
is a measurable space, and /i is a a-finite and non-atomic measure. In that case, 
([3]) implies that 

DtF = nin-i if (•, t)) , DtG = m/™_i {g {-,1)) 

and we have 

{DF,DG)l^nm f In-i{f {-.t)) Fn-i{g {-^t)) ii{dt) . 

J A 

Thanks to the multiplication formula for multiple stochastic integrals, see, for 
instance. Proposition 1.1.3. in 6], one obtains 

/nAm-l / _ 1 \ / _ 1 \ 
E r\\^ ^ ^ jIn+m-2-2r{f{;t)^rg{-,t))fi{dt). 

Taking into account the orthogonality between multiple stochastic integrals of 
different order, wc have 



L[{DF,DG)l 

iAm—1 

E ( 



nAm— 1 / -,\2/ -i\2 

n— 1\ /TO— 1 



X / {f{-,t)§)rg{-,t),f{-,s)§)rg{-,s))He(r.+m-2-2.)fi{dt)fi{ds). 
Notice that 

/ f{-,t)®rg{-,t)^ji{dt)^f®r+i9. 

J A 

and, as a consequence, 

f {■,t)®r9{-,t) ^{dt) = f®r+ig- 

J A 

Therefore, 

nAm— 1 , ,\2/ 

n{DF,DG)l]=n^m^ J2 i^'-fy,) ||/^^+i5||io(„-.™-.-.., , 

which implies the desired result. ■ 
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3 Main result 

Fix n > 2, n e N. Consider a sequence of square integrable random 

variables belonging to the nth Wiener chaos. We know that 

E[\\DFk\\l]=n\\f4l^„. (5) 



The next lemma establish the equivalence between the convergence in (fl) of 

\\DFk\\l toac 
same constant. 



||iDJfc||^ to a constant, and the convergence of E[||i)i^fe||^] to the square of the 



Lemma 3 Consider a sequence {Fk = /„ (/fc)}fcgN of square integrable random 
variables belonging to the nth Wiener chaos such that 

E[F^]=\\Mle.^ 1. 

The following conditions are equivalent: 

1. \iuik^+^E[\\DFu\\H]^n^- 

2. limfe^+ec^ llmll^ - n, m (17) . 
Proof. Notice that, using ([5]), 

m\DFk\\l-nf] = n\DFu\\H]-2nn\DFk\\\]+n'' 

4 1 0„2 II r ||2 , 2 



n 



and the result follows easily. ■ 

Now, we establish the main result of this paper. 

Theorem 4 Consider a sequence {F^ = In of square integrable ran- 

dom variables belonging to the nth Wiener chaos such that 

^[F'^\^\\fk\\le. 1. (6) 

/C— »+CXD 

The following statements are equivalent. 

i) As k goes to infinity, the sequence {-Ffe}j,gN converges in distribution to the 

normal law (0, 1). 

ii) limfc^+eoE[F4] = 3. 

iii) For all 1 < I < n - 1, limfe^+oo \\fk ®i fk\\H»^(^-n = 0. 

iv) \\DF,\\l f-^'' n. 
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Proof. We will prove the following implications 

iv) ^ i) ^ ii) => iii) => iv). 

[iv) ^ i)] The sequence of random variables {-Ffcl^gN is tight because it is 
bounded in L^(f2) by condition Then, by Prokhorov's Theorem we have 
that is relatively compact, and it suffices to show that the limit of any 

subsequence converging in distribution is A^(0, 1). Suppose that, for a subse- 
quence {fcijjgN C {fcjj^^gpj we have 

Fk, , ^ G. (7) 

By condition (O G G (il). Therefore, the characteristic function Lpit) = 
£[6***^] is differentiable and (p' {t) = iE[Ge^*-^]. For every fe G N, define tpk {t) = 

By the Continuous Mapping Theorem, ([7]) implies that 

Ffc^e"^*-- ^ Ge**^. (8) 

I — ^+oo 

The boundedness in L'^iQ) plus the convergence in law ([5]) imply convergence 
of the expectations. Hence, we obtain 

V'k, it) , -> ^'(t). 

On the other hand, using the definition of the operator L, Proposition [1] and 
the definition of the operator 5, we have 

= ^E[{DFk,D{e''^''))H] = ^E[e"^'= \\DFk\\l]. 



Therefore, 



Furthermore, 

Eie^'^-i \\DF,Ml] - n^{t)\ < E[| |1^ - n\] + n |E[e'*^'=.] -^{t)\ 

which, by the definition of ip and hypothesis iv), gives that 

^'k, it) , ^ -tip it) . 

This implies that ip (t) satisfies the following differential equation 

p'{t) = -tpit) 
piO) = 1, 
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which is the differential equation satisfied by the characteristic function of the 
iV(0,l). 

[i) ^ ii)] and [ii) ^ iii)] These implications are proved by Nualart and Pec- 
cati in Proposition 3, [J. The proof of the first one is trivial and the proof of the 
second one involves some combinatorics and it is based on the product formula 
for multiple stochastic integrals. 

[iii) iv)] By Lemma|3]it is enough to prove that iii) implies limfc^+oo 
n?. Using ^ we obtain 

It follows that E[||Z?Ffe||^] converges to if and only if 

||/fe<i'i/fc||!,02(„-i) 0, l</<n-l. 

fc— f-+oo 

As 

||/fe^;//c||^02(„-i) = (2(rt - /))! ||/fc«';//c||^(S2(„-i) < (2(rt - 0)! IIA- /fcll^cs2(.-i) , 

we conclude the proof. ■ 

As a consequence, we obtain. 

Corollary 5 Fix n > 2 and F belonging to the nth Wiener chaos such that 
K[F^] = 1. Then the distribution of F cannot be normal andE,[\\DF\\^] ^ . 

Proof. If F had a normal distribution or E[||_DF||^] = n^, then, according to 
Theorem H] and Lemma [3l we would have Var[||Z?F||^] = 0, but this implies 
= or belonging to the first chaos. ■ 

4 Multidimensional case 

In this section we give a multidimensional version of Theorem 21 For d > 2, fix 
d natural numbers 1 < Hi < ■ ■ ■ < n^. Consider a sequence of random vectors 
of the form 

= {Fl, ...,Ft) - [fl) , ...,/„, {ft)) , (9) 
where fl € iy0"« ^ and 

Ik = ilk')l<^.J<d = i{DFl,DFl)H)l<^.J<d. 

The following lemma shows that the convergence of the covariance matrix of 
Fk to a diagonal matrix plus the convergence in (fl) of the diagonal elements 
of 7fe to the constant n.^, implies the convergence in (f2) of 7^ to a diagonal 
matrix. 
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Lemma 6 Let {-Ffej^gpj he a sequence of random vectors as 0) such that, fo 
every I < i, j < d, 

where Sij is the Kronecker symbol. We have that 



lim E[^^J^]=<5, 

k — ^+oo 



implies 



I ''llfffe^+oo - - 



Ik ^-^^ y/niUjdij,! <ij < d. 

k — *+oo 



Proof. We need to show that, for i < j, one has 



lim E[{DFlDFl)l] = 0. 



Using Q we obtain 

71 i 

E[{DFl,DF^r„] = 



^(K-r)!K--r)!(r-l)!r 



(10) 



(11) 



< 



(njlrij!) (ui + Uj — 2r)! 
^(K-r)!K--r)!(r-l)!f 

We have reduced the problem to show that 



fk^rfk 



fk fl 



2 



lim 

k — >+oc> 



fk^rfl 



H 



»(»ii+n,--2r) 



= 0, 1 < r < n,. 



The next step is to relate the norm of ®r fl with the norms of fl < 
and fl (X)„j--r Using the definition of the contractions, we have 



yyii — r J 



/fc ®r fl ^«(„^ + ,.^_,., - (/fc ®n^~r fl fl ®n,-r fl) H^i-^> 



Hence, by Cauchy-Schwarz's inequality, we obtain 

.2 

fk '^r fl ^^(„.+„_2r) - '^rii-r /fc||^«(2r) 



k '^rij-r 



fk 



(12) 



(13) 



In the case I < r < Hi — 1, hy assumption (jlip and Theorem [4] (implication 
iv) =5> iii)), the right hand side of equation (fT3|) tends to zero as fc tends to 
infinity. In the case r = < Uj, we have that the right hand side of equation 
(fT51) is equal to 



k "yuj-r Jk 



fk 



i/®(2r) 



which tends to zero as k tends to infinity, because 

II ' 1 1 2 
sup /fe L»„, < +00, 
k>l 
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thanks to assumption (jlOp and, analogously to the previous case, 



In the case r = n,- 



the equahty (fTO]) gives 



fl ®r fl 



H' 



3("i + ",--2r) 



which tends to zero by assumption ([TT 

Let Vd be the set of all (ii, 12, is, «4) G (1, ...,d) , such that one of the fol- 
lowing conditions is satisfied: (a) ii ^ 12 = is = *4, (b) ii 7^ 12 = *3 7^ H and 
*4 7^ *i, (c) the elements of (ii, «2, is, ^4) are all distinct. 

The following is a multidimensional version of Theorem [H 

Theorem 7 Let {^fc}/jgpj &e a sequence of random vectors of the form such 
that, for every 1 < i, j < d, 



hm E[FlFl]=S,„ 

k — *+oo 



(14) 



where 6ij is the Kronecker symbol. The following statements are equivalent. 

i) For every j — l,...,d, F^, converges in distribution to a standard Gaussian 

variable. 

ii) For every i = 1, ...,d, limfc^+oo E[(F^)''] = 3. 

iii) For alll<i<d,l<l<n,- 1, \\fl /^IP «,(„^_,) ^ 0. 

iv) For alll<i<d, \\DFl\W ^^^^^ n,. 

v) For every (ii, 12, is, ^4) e Vd, 



lim E 

k — )-+oo 



and 



hm n\{Fl! 



k — >-\-oo 



0. 



vi) As k goes to infinity the sequence {-FfejfcgN converges in distribution to a 
d-dimensional standard Gaussian vector Nd {Q,Id)- 
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Proof. As in the one-dimensional case, we provide a proof to the above theorem 
using Malliavin calculus and avoiding the Dambis-Dubins-Schwarz theorem. 
The equivalences i)<S=^ ii) <S=^ iii) iv) follow from Theorem [D The fact 
that vi)=4' v) is easy, and the implication v)^=> iii) is proved by Peccati and 
Tudor in \10\ using the product formula for multiple stochastic integrals. Hence, 
it only remains to show the implication iv)=^ vi). 

The sequence of random variables {-Ffe}i,gj^ is tight by condition p4|) . Then, 
it suffices to show that the limit in distribution of any converging subsequence 
{Ffcjjgj^ is NdiOJd)- For every k G N, define ipkit) = E[e*<*'-^'=>]. Let ip{t) be 
the limit of (pki (t) as I tends to infinity. As in the proof of Theorem 2] we have 
for all j = 1 , . . . , d 

^(t) ^ ^(t) 

On the other hand, using the definition of the operator L, Proposition [1] and 
the definition of the operator S, we have 



■> h=l 



Therefore, 



^'j^{t) = -^^Y.'nne'^'-'^^h'::]- (15) 

^ "-^ h=l 



Using Lemma [5] and taking the limit of the right-hand side of expression (|15p 
yields 

for j = 1, . . . , d. As a consequence, ip is the characteristic function of the law 



5 Central limit theorem for square integrable 
random variables 

In this section, we will establish a weak convergence result for an arbitrary 
sequence of centered square integrable random variables. 

Theorem 8 Let {-Ffe}j,gpf be a sequence of centered square integrable random 
variables with the following Wiener chaos expansions 

oo 
n=l 

Suppose that 
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(i) liniAr^+oo limsupfe^+^ J2'^=n+i E[(J„-Ffc)^] = 0, 

(ii) for every n>l, limfe^+oo E[( J„Ffc)2] = al, 

(iii) J2n=i = < +00. 

(iv) for all n> 1, 

II n / 7 IP M|2 ^'(") 2 
\\D{JnFk)\\H , n<j„- 

Then, converges in distribution to the Normal law N{0,a'^) as k tends to 
infinity. 

Proof. By Theorem 21 conditions (ii) and (iv) imply that for each fixed n > 1 
the sequence { Jri-Ffc}j,gpj converges in distribution to the normal law N (O, cr^) . 
Moreover, by Theorem [7] we have, for each n > 1, the following convergence 

{JiFk,...,JkFk) ^ ^ (6, -,?«), (16) 

where {Cn}„gN independent centered Gaussian random variables with vari- 
ances {<^n}n^fi- Foi' every > 1, set 

N 



Fk = X! ■ 

N 



Define also ^ = X^^ifn- Let / be a function such that |/| and |/'| are 
bounded by one. Then 



< 



< 



\nf{Fk)]-nfm\ 

|E[/(^^,)] -E[/ (Ff )]| + |E[/ (Ff )] -E[/(e^)]| + |E[/(e^)] -E[/(0]| 

/ 00 \ 1/2 

E E[(J„Ffe)M +|E[/(Ff)]-E[/(e^)]| + |E[/(e^)]-E[/(e)]| 

\n=Af+l / 



Taking first the limit as k tends to infinity, and then the limit as N tends to 
infinity, and applying conditions (i), (iii) and [16] we finish the proof. ■ 



Remark 9 If Fk G ID)^'^ and 



supE(pFfc|l2,) <oo, 

k 



then condition (i) holds. 
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This theorem can be apphed to random variables not belonging to D"'^'^, 
and it requires the convergence of the derivatives of the projections. In this 
sense, it would be interesting to study the relation between the convergence in 
distribution of a sequence Fk to a normal law and the convergence in L^(rj) of 
ll-^-^fell// t;o a constant. As we have seen, these conditions are equivalent (if F). 
are centered and with limfe^oo = a^) for random variables in a fixed 

chaos. In the general case, we conjecture that this equivalence does not hold. 

6 Example 

Suppose that B" = {B^^ ,t > 0} is a fractional Brownian motion (fBm) with 
Hurst parameter H G (0, 1). That is, B^ is a Gaussian stochastic process with 
zero mean and the covariance function 

E(i?f ) = 1 (^2^ + s^" -\t- s\^") . 

Fix H < ^ and an odd integer k > 1. We are interested in the asymptotic 
behavior of 

[nL] 

as n tends to infinite, where t e [0, T]. Set Xj = Bf - Bf_^. Then, {Xj,j > 1} 
is a stationary Gaussian sequence with zero mean, unit variance and correlation 
pH{n) — H(2H — l)n^^~^ as n tends to infinity. We have the following result. 

Theorem 10 The two-dimensional process (B^ , Z^"'^) converges in distribu- 
tion in the Skorohod space X'([0, T])^ to {B^,cW), where W is a Brownian 
motion independent of B^ , and 

oo 

c' = Y.^[{x,x,+jr]. 

Proof. The proof will be done in two steps. 

Step 1. We will first show the convergence of finite-dimensional distri- 
butions. Let {ak,bk\^ k = 1,...,N, be pairwise disjoint intervals conatined 
in [0,T]. Define the random vectors B = (B^ - B^^,...,B^^ - B^J and 

X(") = (x("\...,x(r)), where 

[nak]<j<[nbk] 

We claim that converges in law to {B,V), where B and V are inde- 

pendent and y is a Gaussian random vector with zero mean and independent 
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components with variance c^ibk — Ofc)- By the self-similarity of the fBm, it 
suffices to show the convergence in distribution of F^")) to {B, V), where 



B 



in) 



in) 



[na^]<j<[nbk] 

7^ ^ 

lnak]<j<lnbk] 



with Xj = Bf - Bf_j^ and 1 < k < N. 

We denote by H the Hilbert space defined as the completion of the step 
functions on [0, T] under the scalar product 

(l[o,t]- l[o,s])^ = ). 
Using Theorem [8] it suffies to show that: 



hm E(b(")b["^ 
hm Elsf'j^yf) 

hm E(j„,r,(")j„r,(") 

n^oo \ 



^ E(K-<,)K-BfJ), (17) 
= 0, (18) 
= 4/i(fefe - afc)c^m' (19) 



and 



lim 



(„) 



H 



{bk ~ ak)m(T 



2 

in ' 



(20) 



in L^, for all 1 < to < k, and 1 < k,h < N. The variances cr^ must satisfy 
J2m=i = c^- The projection on the mth Wiener chaos JmYk"^ has the form 
J2[nak]<j<[nbk] Hm{Xj), where Hm{x) denotes the mth Hermine polynomial. 
The convergences (|17l) and (IT51) are immediate. To prove ([T5)) we write 



J2 ^{Hra{X^)Hm{X,)) 

["■ifcl<3<["i>fcl 

OO 

SkhC^nih - flfe) V —:PhU), 



as n tends to infinity. On the other hand, we have 



(n) 



n 



[nafc]<j<[n6fc] 



^ V i/„,_i(XOi/m-i(X,)pH(j-i) 
n ^ — ' 

lnak]<ij<lnbk] 
2 [nbk] [n6fc]-[nafc] 
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^ ^ H^_i{X,)H^_i{X,+,)pH{j). 



i=[nak] + l j=0 
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We claim that the series 



EOO TT 
7=0 -f^" 



-i{X.i)H„i^i{Xi+j)pH{j), converges al- 



most surely and in L^, and {^i, i > 1} is a stationary ergodic sequence. The con- 



L {Xi)Hm- 1 {Xi+j ) I ^ 



< 



vergence in follows from the fact that sup^ E 

oo, and J2'jLo \PH{j)\ < oo. On the other hand, the sequence > 1} is 
ergodic becauuse {Xi,i > 1} is so. Hence, by the ergodic theorem we have in 
L2 



lim 

n— >cxD 



(") 



n 



a-k 



OO 

y 



1) 



■phUY 



which implies (^0)1 . 

Step 2. Taking into account that all norms, for 1 < p < oo, are equivalent 



on a fixed sum of Wiener chaos, in order to show that the sequence 
in I?([0,T])Mt suffices to show that 



(n) 



is tight 



E( 



Z. 



(n) 



-E 



[nt] 

E 

= [ns] + l 



< c\t- 



and this follows easily as above. ■ 

The convergence of the finite dimensional distributions in the above theorem 
can also be deduced from general central limit theorems for functionals of Gaus- 
sian stationary sequences satisfying the Hermite rank condition (see Brauer and 
Major [1]). A related result for the function g{x) = \xY' — E(|i?f^|''), where 
p > and e (0, 1) was obtained by Corcuera, Nualart and Woerner in [5]. 
The central limit theorem was proved in this case using the approach of Nualart 
and Pcccati [7] (see [2], Proposition 10). 

The above result is motivated by the extension of the Ito formula to the 
fractional Brownian motion in the critical case H = ^, using discrete Riemann 
sums. This problem has been considered by Swanson in [11 in the case of the 
solution of the one-dimensional stochastic heat equation driven by a space-time 
white noise {u{t, x),t > 0}, which behaves as a fractional Brownian motion with 
Hurst parameter H = j, for any flexed x G M. In this case, the convergence 
in law to a Brownian motion is proved for a modiflcd sum of squares of the 
increments. 

Acknowledgement 11 David Nualart would like to thank Jason Swanson for 
very stimulating discussions regarding the Example in Section 6. 
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